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The renormalization group method is applied for obtaining 
the asymptotic form of the wave function of the quantum 
anharmonic oscillator by resumming the perturbation series. 
It is shown that the resumed series is the cumulant of the 
naive perturbation series. Working out up to the sixth order 
and performing a further resummation proposed by Bender 
and Bettencourt, we find that the agreement with the WKB 
result becomes worse in the higher orders than the fourth at 
which the agreement is the best. 

PACS numbers: 03.65.Ge,02.30.Mv,H.15.Bt,11.15.Tk 



It is well known that naive perturbation series are di- 
vergent or at best asymptotic series (|] . One needs to re- 
sum the divergent series to obtain a sensible result from 
the perturbation theory. Indeed, various resummation 
techniques have been devised [[[) . Recently, a unified and 
mechanical method for global and asymptotic analysis 
has been proposed by Goldenfeld et al j| . This is called 
the renormalization group (RG) method. The unique 
feature of their method is to start with the naive pertur- 
bation theory and allow secular terms to appear in con- 
trast with all previous methods Q; adding unperturbed 
solutions to the perturbed solutions so that the secular 
terms vanish at a "renormalization point" t = to and 
then applying the RG equation, one obtains a resummcd 
perturbation series. 

Subsequently, the present author formulated the RG 
method geometrically on the basis of the classical theory 
of envelopes || ; it was indicated that the RG equation a 
la Gell-Mann-Low in the field theory can be identified as 
the envelope equation. 

The purpose of the present work is to apply the RG 
method as formulated in j3j to Schrodinger equation of 
the quantum anharmonic oscillator (AHO) and obtain 
the asymptotic form of the wave function. 

AHO is a theoretical laboratory for examining the 
validity of various approximation techniques 0,^]. Re- 
cently, Bender and Bettencourt || have shown that 
multiple-scale perturbation theory (MSPT) can be suc- 
cessfully applied to the quantum anharmonic oscillator; 
MSPT or the reductive perturbation theory is known to 
be one of the most general method in applied mathemat- 
ics 0| apart from the RG method for improving pertur- 
bative expansions |?]] . They examined the Heisenberg op- 
erator equation and the (time-independent) Schrodinger 
equation: The exact closed-form solution was found for 



the Heisenberg equation, and the asymptotic behavior 
of the wave function ip(x) for large x was constructed 
which agrees with the WKB result. One should remark 
here that a further resummation had to be adopted for 
the latter case, which is not intrinsic in MSPT and a simi- 
lar method had been proposed by Ginsburg and Montroll 

Actually the RG method a la Goldenfeld ct al has 
been already applied to quantum mechanics by others 
H|o|]: Egusquiza and Valle Basagoiti (|] applied it to 
solve the time-dependent Heisenberg equation considered 
in ||, while Frasca used it to solve the time-dependent 
Schrodinger equation of a two-level system. There has 
been, however, no attempt to apply the RG method to 
obtain asymptotic forms of wave functions. The rea- 
son may be that it is not trivial to identify the secular 
terms for the wave functions which can be made to van- 
ish at a "renormalization point" x — xq. In the present 
work, such secular terms are successfully identified for 
the ground and the first excited states. 

We shall show that the resummation of the perturba- 
tion series of the wave functions is performed in the RG 
method more mechanically and explicitly than in MSPT: 
The polynomials f n (x) in the resumed series are explic- 
itly given in terms of the polynomials P n (x) obtained in 
the naive perturbation theory Furthermore, it will 
be found that ^2^Lq £ n fn(x) is the cumulant series pd| ] 
of the naive perturbation series J2^=o eV "Pn(x)- Since our 
method is mechanical and easy to perform, we shall work 
out up to the sixth order in the perturbative expansion 
and examine how the results of Bender and Bettencourt 
persist or are modified in the higher orders. 

Our Hamiltonian for the anharmonic oscillator is given 

by @ 



H = p 



1 2,1 4 

-x + -ex , 



4 4 

and we consider the following Schrodinger equation 
(H-E)tP(x)=0, 



(1) 



(2) 



with the boundary condition ?/>(±oo) = 0. We shall con- 
fine ourselves to the ground state for the moment. WKB 
analysis shows that for large x, 



ip(x) ~ exp{-v^|a;| 3 /6}. 



(3) 



We shall examine how the perturbation theory can re- 
produce the WKB behavior or not, as was done in |6j. 
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As preliminaries, we first apply the Bender- Wu method 
[ H for performing Rayleigh-Schrodinger(RS) perturba- 
tive expansion. The wave function and the eigenvalue are 
both expanded as power series of e; 

oo oo 

^(i) ~ £Vy n (o;) and E{e) ~ £ e n E n . (4) 

n=0 n=0 

We take the boundary values hence the normalization as; 

Ito(O) = 1 and y„>i(0) = 0. (5) 

The lowest order solution reads 

yo (x) = e -* 2 /4 an d B = 1/2. (6) 

The higher order terms with n > 1 are written as 

y n (x) = er'^^P^x), (7) 

where P n (x) is a polynomial. It is readily shown that the 
polynomials satisfy the recursion relation; 



n-l 



J*'(x) - asi*; (z) - T P„-i(x) - J2 p 3 ix)E n -j ■ (8) 



2;; 



fc=i 



where the coefficient C n k satisfies a recursion relation 



2fcC„, fc + C n _i >fc _ 2 = -(A + l)(2fc + 1)C, 



n-l 

£ 



n,fc+l 



Cj.feCn-j',!) (11) 



Now we apply the renormalization group method to 
resum the perturbation series obtained above. We shall 
present the method so that it becomes clear that the no- 
tion of envelopes is intrinsically related to the method 
H : We shall also make it clear that the RG method con- 
cerns with the boundary conditions in conformity with 
the general property of the RG methods as emphasized 
by Shirkov @. 

First we try to obtain the wave function ip(x;xo) 
around an initial point x — xq in a perturbative way: 



ip(x; xq) ~ e n z n (x; xq) and E(e) 



e n E n , (12) 



with the initial or boundary condition (BC) at x — xq; 



tp(x ;x ) = W(xq). 



(13) 



We suppose that the boundary value W(xo) is always on 
an exact solution of Eq.(2). W(xo) may be also expanded 
in a power series of e 



W(x ) 



l W n (x ) 



(14) 



This equation determines the polynomials and the eigen- 
values E n successively. Here it should be remarked that 
y n {x) in Eq.(7) may be identified as a secular term be- 
cause it is a product of the unperturbed solution and a 
function that increases as x goes large. 

The eigenvalues E n are given by the condition (solv- 
ability condition) dx yo(x)hF(x) — 0, where h = 
(d/dx) 2 — xd/dx and F(x) is an arbitrary function with 
which the integral converges. Thus one finds 

1 f°° x 4 ^ 

E n = 7n= \ dx yo(x){-rP n -i(x) - V P j (x)E n _ j }. 

2 V 7T J-oo 4 ~{ 

Note that E n is determined in terms of only the polyno- 
mials Pj(x) with j < n — 1. 

With these eigenvalues Ej (j < n), the polynomial 
P n [x) is determined by Eq.(8). The general form of the 
polynomial is expressed as 



If we stop at, say, the iV-th order, we will have 
J2n=o e7lz n(x) = ^ N ^(x; xo) which is valid only locally 
at x ~ xq. However, one may take another point of view 
as follows: Geometrically, we have a family of curves 
{i/j( n \x;xo)} Xo parametrized with xq, and each curve 
of the family is a good approximation around x = Xq. 
Then, if each curve is continued smoothly, the resultant 
curve will valid in a global domain of x. This is noth- 
ing else than to construct the envelope of the family of 
curves. More specifically, we only have to determine the 
/gxboundary values W n (xo) so that the perturbative solu- 
tions around x = xo form an envelope. This is the basic 
strategy of the RG method described geometrically. Fur- 
thermore, to be as accurate as possible, the lowest value 
Wq(xq) should approximate the exact value ip(xo,Xo) as 
close as possible, or W n >i(xo) should be made as small 
as possible. 

Let us perform the above program. First, we note that 
the lowest order solution may be written as 



(10) 



z (x;x ) = A(x Q )e l2/4 , 



(15) 



and Eq = 1/2; we have made it explicit that the ampli- 
tude A(xq) may be dependent on Xq. The choice of the 
lowest order solution implies that we have also chosen the 
boundary value as 



Wo(x ) = A(x )e- x o/\ 



(16) 



with C n> i — E n . The recursion relation is solved for low 
k with given n. The polynomials P n (x) up to the six 
order are presented in H, which we refer to. 



The higher order terms with n > 1 may be written as 

z n (x; x a ) = A(x )e~ x2/4 Q n (x] x ), (17) 
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where Q„(x;xo) is a polynomial of x, dependent on xq. 
It is readily shown that the polynomials satisfy the same 
recursion relation Eq.(8) as P„(x). However, since we 
want to make the boundary value Wo(xq) as close to the 
exact one as possible, we impose the boundary condition 
as 

z„(x ;x ) = W n (x ) = or Q n (x ; x ) = 0, (18) 
for n > 1. 

Since Qi(x;x ) satisfies the same equation as Pi(x) 
does, one readily obtains 



Thus we obtain the approximate solution valid around 



x <~ Xq 



tp(x;x ) ~ A(x )e * 2/A ^ e"Q„(x; x ), 



n=0 



which satisfies the boundary condition 

ip(x ;x ) = W Q (x ) = A(x Q )e' x o /4 . 



(25) 



(26) 



l(x;x ) = Pi(x) - Pi{x ), 



(19) 



which satisfies the BC Eq.(18). Notice that a constant is 
the solution of the homogeneous equation. The second 
order equation now reads 

x 4 2 

Q 2 (x;x ) - xQ 2 (x;x Q ) = (—P 1 (x) - ^P j {x)E n ^ j ) 

3=0 

- Pl ( Xo )(^- El ). (20) 

One can verify that E 2 is given Eq.(9). Since Eq.(20) is 
linear and the inhomogeneous part is a linear combina- 
tion of those for P 2 (x) and Pi(x), Q2(x;xo) is given by 
a linear combination of Pi (x) and P\ (x) ; 

Q 2 (x;x ) = (P 2 (x) - P 2 (x )) - Pi(x )(Pi(z) - Pi(x )), 

which satisfies the BC Q2{xq;xq) = 0. One finds that 
Qs(x;xo) satisfies 

x 4 3 

Q'l{x;x Q ) -xQ' 3 (x;x ) = {-^(a:) -^Pji^n-j} 

3=0 

4 2 
-P 1 (x ){^P 1 (x)-J2P3(x)En- J } 
3=0 

-(P 2 (.*o)-Pi(* ) 2 )4-i?i). (22) 

One can see that the inhomogeneous part is again com- 
posed of a linear combination of those for P n (x) (n = 
1,2,3). Thus the solution satisfying the BC is found to 
be 

Q 3 (x;x ) = P 3 (x) - P 3 (x ) - Pi{xq)(P 2 (x) - P 2 {x )) 

-(P2(x ) - Pi(x ) 2 )(Pi(x) - Pi (so)), (23) 

We remark that E% is the same as what obtained for 
P 3 (x). 

Repeating the procedure, one finds that Q n (x;xo) arc 
expressed in terms of Pj (x) (j < n) . For instance, 

Q 4 (x;x ) = P 4 {x) - P 4 (x ) - Pi(x )(P 3 (x) - P 3 (x )) 
-(P 2 (zo) - PiK) 2 )(P 2 (x) - P 2 (x )) 
-(P 3 (x ) - 2P 1 (x Q )P 2 (x Q ) + Pi(x ) 3 ) 
■(Pi^-Pr^o)). (24) 



One may say that now we have obtained a family of 
curves {ip(x; Xq)} Xo with xo parametrizing the curves. If 
x < x < x + Ax with Ax being sufficiently small, the 
wave functions ip(x;x n ) and tp(x;x n + Ax) should give 
the same value at x, i.e., 



ip(x; x ) = ip(x; x + Ax). 



(27) 



Taking the limit Ax — > this condition is found to yield 
that 



dip(x;x ) 



dxo 



= 0. 



(28) 



Notice that when Ax — > 0, x — > xq. This is the basic 
equation of our method. This is nothing but the condi- 
tion to construct the envelope of the perturbative wave 
functions valid around x ~ xo. It is apparent that the 
(21e)quation has the same form as the renormalization group 
equation, hence the name of the RG method. The equa- 
tion gives a condition which A(x ) must satisfy; 



"' 4 =A(x)-j-- f>"(-Q„(x;x )) 
dxo ' 



dx 



n=0 



Defining f n (x) by 



d r, i n| df n {x) 



one obtains 



A(x) = A-exp{J2 e n fn{x)]. 



(29) 



(30) 



(31) 



With this solution, the global solution iI>e(x) is given by 
the boundary value by construction; 



i> E (x) = W Q (x) = A(x)e- 



2 /4 



= A- + f n (x)]. (32) 

71=1 

This is one of the main results of the present paper. 

/„(x)'s are easily calculated in terms of P n (x), and we 
have 



3 



fi(x) = Pi(x) = ~x 2 -±x\ 

Mx) = P 2 (x) - \pi{x) 2 = §x 2 + ^ + i/, 

f 3 (x) = P 3 (x) - P x {x)P 2 {x) + iPiW 3 , 



333 2 45 4 
-x - —x 



32 



32 



— X 6 - —X* 

192 256 
1 



h{x) = P 4 (x) - P 1 {x)P 2 (x) - ^(x) 2 + P 1 {x) 2 P 2 {x) 



on the construction of an envelope certainly resumes the 
perturbation series of the wave function and the resultant 
expression are given in terms of the cumulants of the 
naive perturbation series. 

Now let us examine how the WKB result Eq. (3) can be 
constructed from the perturbation series obtained above. 
Bender and Bettencourt found that if all terms beyond 
1/512 • e 4 x 10 are neglected, the sum of the highest power 
terms in fj(x) (j < 4) is nicely rewritten as 



30885 2 8669 4 1159 6 



— (1 + 2ex 2 + — e^ 4 + — e 6 x 
4 v 12 12 



1152 



256 



512 



768 



163 
2048 ^ 



„io 



512' 



which behaves for large x as 



fe{x) = P 5 (x) - Pi(x)P 4 (x) - P 2 (x)P 3 {x) + P x {x)P 2 {xY 

-P 1 (x) 3 P 2 (x) + P^xfP^x) + iPiOr) 5 , 

5 

916731 2 33171 , 6453 , 
— x l ~r~rz. — x rrr-^ar 



512 



128 



256 



823 8 319 in 
-x s a; 10 



512 



5120 



6144 



f 6 (x) = P 6 (x) - P x {x)P h {x) - P 2 (x)P 4 (x) + P x {x) 2 P A {x) 
~P 3 (x) 2 + 2P 1 {x)P 2 (x)P 3 (x) - Pi(x) 3 P 3 (x) 

+ l P ^ 3 - l P l(x) 2 P2(x) 2 + Pl{xfP 2 {x) 

65518401 2 19425763 4 752825 fi 
-x 2 H — — — x 4 + x b 



- v^M 3 /4(1152/47) 1/8 ~ Ve|cc| 3 /5.96663 (36) 

in an excellent agreement with the WKB result. How 
about the higher orders. In the fifth order, the sum of 
the highest powers may be rewritten by neglecting all 
terms beyond 7e 5 x 10 /1286 as 



-— • (1 + ex 2 /4 - e 2 x 4 /24 + e 3 x 6 /64 - e 4 x 8 /128 
-7eV°/1286) 



4459 



io\i/io 



2048 
43783 



4096 



1536 



3481^2255^^ (33) 



4096 2048 24576 4096 

and so on. fi(x) ~ f 3 (x) coincide with the results in 0, 
where explicit expressions of f n (x) are given only for n < 
3. It is interesting that the polynomials f n (x) are given 
in terms of P n {x) appearing in the naive perturbative 
expansion in a closed form. 

Here learned readers may have suspected that fi(x)'s 
(i = 1,2,3,...) are the cumulant [ p~T] | of the sum 
J2lT=o en Pn{ x ) m the sense that 



164608 ' 
For large x, the coefficient of — v^M 3 is 
4(164608/4459) 1/10 ~ 5.73827, 



(37) 



(38) 



which deviates from 6 more than the fourth order result. 
The sixth order becomes worse: The sum of the highest 
powers is rewritten as 



4 

67621 



1 + Sex 2 



9 3 6 577 4 8 
8 4 768 



^ 1324349 1/12 



(39) 



n=Q 



e"P„(x)~exp[V e n f n (x)] 



71=0 



(34) 



In fact, this is the case. When a function C(£) of £ is 
given by C(£) = X)^o^™/ n ' ' ^ ne n ~th cumulant 
A„ is defined by lnC(£) = Xm^o^"/ 71 ' ' Ki- Expanding 
lnC(£) = ln(l + £(mi + £ 2 /2 ' M2 + •••). one nnds that 
Ai = hi, X 2 = M2-M?! ^3 = M3-3/iiA*2 + 2juf, A 4 = ju 4 - 
4/^i/i 3 -3^+12^ 2 ^ 2 -6^ 4 , A 5 = 24/xf— 6Q / uf/x 2 +30yLti / u|+ 



20^if^3 -10^2^3-5^1^4 + ^5, A 6 



-120^ + 36 °MI^2- 



270fiit4 +30(4- 120^3 + 120^2^3 -10/x| + 30//f// 4 - 
15/i2M4 — 6^i^5 + /ig, and so on. Putting fin = ^'-Pn and 
A n = n\f n , one sees that the relation Eq. (33) between 
P n and /„ is reproduced. In short, the RG method based 



493824 35555328 
which makes the coefficient of — i/e|a;| 3 for large x 

4(35555328/1324349) 1/12 ~ 5.26181. (40) 

In summary, we have successfully applied the RG 
method as formulated in J3| to Schrodinger equation of 
the quantum anharmonic oscillator (AHO): The naive 
perturbation series of the wave function are resummed 
by the RG equation. We have seen that the resumma- 
tion is performed in the RG method more mechanically 
and explicitly than in MSPT. We have shown that the re- 
summed series X^^Lo e "/n( a; ) is the cumulant of the naive 
perturbation series. We have worked out up to the sixth 
order in the perturbative expansion and found the fol- 
lowing: Although the sum of the highest power in f n (x) 



4 



can be organized so that it becomes asymptotically pro- 
portional to -\/eM 3 as was done in ||||, tne coefficient 
of it reaches the closest value to 6, the WKB result, in 
the fourth order, then goes away monotonously from the 
closest value in the higher orders. This is plausible be- 
cause the convergence radius of the perturbation series is 
zero; the cumulant series should be at best an asymptotic 
series. 

We remark that the RG method as developed here can 
be also applied to the first excited state; it is, however, 
unlikely that the method can be used to the higher ex- 
cited states beyond the first excited state. 

Finally, we mention that a variational perturbation 
method called the delta-expansion method fllj has been 
extended for obtaining wave functions [fl6| . The key in- 
gredient of the extension is to construct an envelope of a 
set of perturbative wave functions as in the RG method 
but with a variational parameter. In this method, al- 
though the basic equation can not be solved analytically 
but only numerically, uniformly valid wave functions with 
correct asymptotic behavior are obtained in the first- 
order perturbation even for strong couplings and for ex- 
cited states. In the present method, the basic equations 
are solved analytically, and the asymptotic form of the 
wave function is constructed explicitly, although a fur- 
ther resummation devised in |^|J^] is needed for obtaining 
the asymptotic form. In this sense, the two methods are 
complementary. 

This work was partially supported by the Grants-in- 
Aid of the Japanese Ministry of Education, Science and 
Culture (No. 09640377). 



Theor. Phys. 97, 1 79 (1997); see also a comprehensive re- 



[2] 



C. M. Bender and S. A. Orszag, Advanced Mathematical 
Methods for Scientists and Engineers (McGraw-Hill, New 
York, 1978). 

N. Goldenfeld, O. Martin and Y. Oono, J. Sci. Comp. 
4,4(1989); N. Goldenfeld, O. Martin, Y. Oono and F. 
Liu, Phys. Rev. Lett. 64, 1361(1990); N. D. Goldenfeld, 
Lectures on Phase Transitions and the Renormalization 
Group (Addison-Wesley, Reading, Mass., 1992); L. Y. 
Chen, N. Goldenfeld, Y. Oono and G. Paquette, Physica 
A 204, 111(1994); G. Paquette, L. Y. Chen, N. Golden- 
feld and Y. Oono, Phys. Rev. Lett. 72,76(1994); L.-Y. 
Chen, N. Goldenfeld and Y. Oono, Phys. Rev. Lett. 73, 
1311 (1994); Phys. Rev. E 54, 376 (1996). 
For recent developments, see for example, R. Gra- 
ham, P hys. Rev. Lett. 76 (1996) 2185; S. Sasa, 
preprint patt-sol/9608008, to be published in Phys - 
ica D; K. Matsuba and K. Nozaki, patt-sol/9702001 , 
p_att-sol/9703004| ; H.J. de Vega and J.F.J Salgado, 



[3] 



tiep-th/9701104, to be published in Phys. Rev. D56. 



[4] 

[5] 

[6] 
[7] 

m 

[9] 

[10] 
[11] 

[12] 
[13] 



view; T. Kunihiro, Datt-sol/9709003, which is a refine- 
ment of the argument given in the previous papers where 
the fact that the RG method is a theory manipulating 
initial conditions were not fully recognized. 
Large Order Behavior of Perturbation Theory, Current 
Physics - Sources and Comments, vol.7, ed. J. C. Le 
Guillou and J. Zinn- Justin, (North-Holland, Amsterdam, 
1990). 

G. A. Arteca, F. M. Fernandez and E. A. Castro, Large 
Order Perturbation Theory and Summation Methods in 
Quantum Mechanics, (Springer- Verlag, Berlin, 1990). 
C. M. Bender and L. M. A. Bettencourt, Phys. Rev. Lett. 
77, 4114 (1996); Phys. Rev. D54, 7710 (1996). 
We are aware that MSPT was applied to quantum me- 
chanical problems earlier by Frasca; M. Frasca, Nuovo 
Cimento, B107,915 (1992); 109,603 (1994). 
C. A. Ginsburg and E. W. Montroll, J. Math. Phys. 19, 
336 (1978). 

I. L. Egusquiza a nd M. A. Valle Basagoiti, Report NO. 



[14] 
[15] 



[16] 



hep-th/9611143 



1996. 



M. Frasca, Phys. Rev. A 56,1548 (1997). 
For example, see R. Kubo, M. Toda and N. Hashit- 
sume, Statistical Physics II Nonequilibrium Statistical 
Mechanics, 2nd ed. (Springer 1995). 

We use the notations of Bender and Bettencourt Q for 
later comparison. 

The fact that the RG in general concerns with initial 
or boundary conditions was emphasized by Shirkov; D. 
V. Shirkov, Russ ian Math. Surveys 49: 5, 155 (1994); 
hep-th/9602024 V. F. Kovalev, V. V. Pustovalov and 



D. V. Shirkov, iep-th/9706056 to be published in J. 
Math. Phys. 

C. M. Bender and T. T. Wu, Phys. Rev. 184, 1231 
(1969); ibid. D7, 1620 (1973). 

The following is some of recent articles from which one 
can trace the enormous number of references related to 
this method; C. Arvanitis, H. F. Jones and C. S. Parker, 
Phys. Rev. D52, 3704 (1995). W. Janke and H. Kleinert, 
Phys. Rev. Lett. 75. 2787 (1995). E. J. Weniger, Phys. 
Rev. Lett. 77, (1996) 2859. H. Kleinert, Path Integrals 
in Quantum Mechanics, Statistical and Polymer Physics, 
2nd. edition (World Scientific, Singapore, 1995). 
T. Hatsuda, T. Kunihiro and T. Tanaka, Phys. Rev. Lett. 
78, 3229 (1997); see also S. K. Kauffman and S. M. Perez, 
J. Phys. A 17, 2027 (1984). 



T.Kunihiro, Prog. Theor. Phys. 94, 503 (1995); (E) 95, 
835 (1996); Jpn. J. Ind. Appl. Math. 14, 51 (1997); Prog. 



5 



